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Possibility for J/"^ suppression in high multiphcity proton-proton colhsions at 

v^ = 7TeV 
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We study J/vP absorption in high muhiplicity proton-proton (pp) colhsions at ^sjvjv = 7TeV. 
We predict a modification of the J/\l/-yield within the UrQMD transport approach, where exphcit 
interactions of the J/'i/ with the surrounding comovers and a prehadronic phase with adjusted cross 
sections and J/\[' melting is included. We present the analog of the nuclear modification factor 
in pp collisions at various charged particle multiplicities. It turns out that J/5'-Mesons may be 
suppressed towards higher particle multiplicities in pp collisions at LHC energies. 
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I. INTRODUCTION 

One major goal of ultra-high energy heavy ion physics 
is to recreate the phase of quarks and gluons (the Quark 
Gluon Plasma, QGP) as it might have existed micro- 
seconds after the Big Bang. Various experimental facili- 
ties have been built to explore the properties of the QGP 
experimentally, while on the theory side a multitude of 
(potential) signatures and properties of the QGP have 
been predicted. Recently, the field of heavy ion physics 
has entered a new era with the start of the Large Hadron 
Colhder (LHC) at CERN. For the first time this offers 
the opportunity to investigate collective effects already in 
proton-proton (pp) collisions, especially when high mul- 
tiplicity interactions are selected [H-Il]- 

In this letter, we explore the absorption of J/^t^-Mesons 
at high charged particle multiplicities in pp collisions at 
^/snn = 7TeV. Similar studies have recently been per- 
formed by I'Jg'] and found a substantial influence of the 
medium at high multiplicities on observables. This study 
is especially interesting, because it i) allows to explore if 
a collectively expanding system has been created in pp, 
ii) can provide insights into the pp baseline necessary 
for the interpretation of heavy ion collisions and iii) pro- 
vides complementary information to the studies of AA 
collisions. Up to now, studies of J/4's in a dense and 
hot medium were restricted to (massive) nucleus-nucleus 
(AA) reactions. Here it was found that the initially pro- 
duced charmonia are suppressed in heavy ion collisions 
compared to pp collisions at the same energy. Three ef- 
fects to explain this suppression have been put forward: 

• Nuclear absorption or baryonic suppression Q. 
These effects depend on the thickness function of 
the two nuclei and therefore change with central- 
ity. This 'normal' suppression with centrality is not 
able to describe the measured J/vP-yield in very 
central AA collisions at SPS and RHIC energies 
Two possible interpretations of this 'addi- 
tional' suppression have been suggested and will be 
discussed next. 



• Comover absorption fl^ - flTj . Here the charmonia 
are additionally absorbed by inelastic scatterings 
with comoving mesons. It is assumed that the cor- 
responding J/^ hadron cross sections are on the 
order of a several mb ilSj . Nevertheless theoretical 
approximations for the absorption cross section dif- 
fer by more than an order of magnitude fl^ . In pp 
collisions at low energies (i.e. low particle densities) 
comover absorption is irrelevant. 

• Debye screening . An alternative explanation of 
the additional suppression is explained by the for- 
mation of a QGP. In a QGP the formation of J/^*- 
Mesons is suppressed due to Debye screening in the 
matter. The charm quarks will leave the reaction 
zone as hadrons with open charm. Therefore J/^ 
suppression in heavy ion collisions has been pro- 
posed as a possible signature for QGP formation 
poj . However, see also [2lj for a recent discussion. 

For massive nuclear collisions at higher energies, it was 
also speculated that J/^'-Mesons might be reformed via 
J/'^f <->■ cc or J/^* -o- DD reactions, if the charm quark 
(or D-Meson) densities become substantial p2| - [2^ . Fur- 
ther information in this respect might also be gained from 
correlation studies as suggested in [2^ [2^ . 

High multiplicity Proton-Proton interactions at LHC 
offer a completely new set of questions and possibilities 
as compared to nucleus-nucleus reactions at similar mul- 
tiplicities. At these energies the particle multiplicities 
at central rapidities rise to the same order of magnitude 
as in heavy ion collisions at lower energies, and the en- 
ergy densities may even exceed the values of central Au- 
Au/Pb-Pb reactions at SPS and RHIC. This allows to 
investigate charmonium suppression at high energy den- 
sities and temperatures in a rather clean environment ex- 
cluding baryonic effects. It does also not suffer from the 
potential charm recombination due to the still low charm 
yields in pp, thus enabling one to explore the absorption 
mechanisms with higher precision. 



Thomas Lang and Marcus Bleicher. Possible J/^ suppression in high multiphcity proton-proton colhsions at LHC energies 2 



II. THE MODEL AND RESULTS 



body transition model [3l|, |32 



For the present study of J/^'-interactions with the 
medium, we modified the UrQMD transport approach 
to include known effects not present in the standard sim- 
ulation. Here we assume a high temperature regime with 
melting of the J/4' (T > 325 MeV), a prehadronic stage 
(155 MeV < T < 325 MeV) and finally a hadronic stage 
(T < 155 MeV). As UrQMD can not directly provide the 
temperature of the medium we map the transition tem- 
peratures to scalar quark densities using the chiral equa- 
tion of state implemented in the UrQMD hydro model 
[27| . We obtain for T 325 MeV a quark density of 
Uq = 12fm~^ and for T = 155 MeV a quark density of 
Uq = 0.6 fm~^. The local scalar quark densities at the 
position of the J/^ are then used to obtain the corre- 
sponding J/^E" cross sections in each UrQMD time step. 
In the following the assumptions for the hadronic, pre- 
hadronic and melting phase are explained: 

Melting regime (T > 325 MeV, > 12/fm^): 
The onset of charmonium melting for J/4's is, on the ba- 
sis of [28J, at 325 MeV. One should keep in mind that 
charmonium states are supposed to still persist in the 
QGP, their wave function should broaden nevertheless 
with increasing energy density (28j . Therefore in the 
present study J/4's are assumed to melt if the energy 
density is high enough and the charmonia stay in this 
hot medium for some time. Here we choose, on the basis 
of [28j, a proper time of 1 fm/c for the charmonium decay 
in the hot medium. 

Prehadronic regime (T > 155 MeV, > 0.6/fm^): 
The temperature of the phase transition is, in-line with 
the lattice QCD estimates [H, [sO], set to 155 MeV. In 
this prehadronic phase we mimic the QGP effects in 
UrQMD by prc-formed hadrons. To calculate charmo- 
nium dissociation we use fixed cross sections in the hot 
medium. The inelastic cross sections to mimic the effec- 
tive dissociation of the QGP on J/^* in the prehadronic 
phase are obtained by adjusting the hadronic cross sec- 
tions to data at SPS energies, see FigH] The reason we 
use SPS data to fix the cross sections is the negligible 
amount of recombination of D-Mesons to charmonium 
states due to the low D-Meson density at SPS energies. 
For this exploratory study we take constant cross sections 
in the prehadronic phase with a meson- J/^" elastic and 
inelastic cross section of 0.78 mb, while the baryon-J/^E' 
cross sections are obtained by quark number scaling. 

Hadronic regime (T < 155 MeV, < 0.6/fm^): 
The hadronic phase includes elastic scattering of J/^- 
particles and its dissociation by baryons and mesons as 
well. Furthermore recombination of D-Mesons to char- 
monium states is implemented using detailed balance, 
but found to be negligible in pp collisions due to the 
low D-Meson abundance. For the inelastic charmonium- 
meson cross sections in the hadronic phase we use a 2- 



fl+2^3+4(s) 



ma -I- 7714 



PL 

Pi 



(1) 

Here Ei denotes the energies of the ingoing and outgoing 
particles. Mi the masses and pi and pf the initial and 
final momenta in the 2-particle rest frame. The effective 
matrix element |Mip — 0.65, is fixed to Pb-Pb colli- 
sions at SPS energies. The matrix element is increased 
by a factor of three for an excited D-Meson in the out- 
going channel and decreased by a factor of three for a 
strange D-Meson in the outgoing channel [3l|, [s^l • The 
corresponding back-reaction for D-Meson recombination 
is given by detailed balance. 



0'3+4^1+2(s) = 0-i+2-i.3+4(s) 



(2^1 + 1) (2^2 



(2^3 + l)(2^4 + l)p^ 



(2) 



where Si denotes the spin of the particles. For the inelas- 
tic baryon-J/vE" interactions we use a constant dissocia- 
tion cross section of 4.18 mb obtained from [^H^. For all 
elastic cross sections of charmonia in the hadronic phase 
we use a constant cross section of 5 mb. In this hadronic 
phase the usual UrQMD formation times are considered 
for the produced hadrons. The formation time of the 
J/^* itself is negligible due to its large mass. 

Let us compare the present parameter set to data ob- 
tained at SPS, which is in the same charged particle re- 
gion as the pp data at LHC. To model the initial state for 
the J/^-particles one needs the momentum distribution 
and the spatial distribution at the production points. 
For the momentum distribution of the produced charmo- 
nia at SPS we use an ansatz from 1341 . 



dN 



dxpdpj 



(1 



\XF 



(3) 



The distribution in a:^ is taken from fsF] with c = a/{l-\- 
hj \fs)^ where a = 16 and h = 24.9 [31|,i36|]. The exponent 
for the transverse distribution is hp^ = 2.08 GeV^^ [33 |. 
Our spatial distribution for the charmonium production 
at SPS is based on the Glauber model and calculated 
using the UrQMD model. Here we perform a pre-run 
for each event where we write down the nucleon colli- 
sion points that we get when we switch off all interac- 
tions. These nucleon collision points are used in the sub- 
following full event for possible charmonium production 
coordinates. Fig. [1] depicts J/ 4' decays to muons scaled 
by the number of Drell-Yan pairs dependent on the num- 
ber of participants Npart and the corresponding number 
of charged particles Nch- UrQMD data are compared to 
data of the NA50 experiment [8]. In line with the ex- 
pectations we observe a strong decrease of the J / yield 
towards central interactions. One should note that the 
present set-up is slightly different from the one suggested 
by Spieles et al. [1^ Eal ■ Here we assume the presence of 
prehadronic states in comparison to the previous analy- 
sis that incorporated a time dependent formation time. 
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FIG. 1: (Color online) J/^ decays to muons scaled by 
the number of Drell-Yan pairs dependent on the num- 
ber of participants Npart and the corresponding number 
of charged particles Nch- UrQMD data are compared to 
data of the NA50 experiment The normalization fac- 

tor B^^cr(J/^')/(T(Dy^.9_4.5 =34.9 in pp interactions at 
158 GeV is taken from 



rapidity 


A 


B 


C 


z 


\y\ < 1-2 


0.459 


6.72 


1.8 


5.04 


1.2 < |y| < 1.6 


0.601 


9.13 


1.8 


8.03 


1.6 < \y\ < 2.4 


0.33 


4.57 


1.8 


4.52 


2.4 < y < 2.5 


0.685 


22.0 


1.1 


21.1 


2.5 < \y\ < 3.0 


0.234 


5.13 


1.7 


5.18 


3.0 < \y\ < 3.5 


0.188 


4.42 


1.8 


4.44 


3.5 < \y\ < 4.0 


0.215 


4.54 


1.8 


4.68 


\y\ > 4.0 


0.282 


4.95 


1.75 


5.46 



TABLE I: Coefficients for the momentum distribution in For- 
mula 5. 
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There are also further differences in the energy depen- 
dence of the hadronic cross sections, which were assumed 
to be constant in the Spieles analysis. 

Let us now turn to the analysis of the LHC pp data. 
The momentum distribution of J/vPs at LHC has been 
measured with the various experiments at a center of 
mass energy of 7TeV. Because we are interested in the 
late stage effects, we do not attempt to model the ini- 
tial state, but fit the rapidity distribution of [37] using a 
double Gaussian fit 



da 



J/* 



dy 



A ( „i 
In ■ a ^ 



(^) 



and the transverse momentum distribution of [36 
[39| using 



dV(J/^) 

dprdy 



Pt 



A- 1 



(4) 



and 



(5) 



with the parameters A,B,C and z. We fit the parame- 
ters for the transverse momentum distribution for each 
measured rapidity bin separately. The parameters are 
displayed in Table HI 

To obtain the initial spatial distribution of the J/^* 
emission points in pp reactions we smear the production 
points located in the proton-proton collision plane with 
a Gaussian distribution of variance a = 0.88 fm. Let us 
now explore pp collisions at a center of mass energy of 
^/sWn = 7TeV within this approach. Since PYTHIA 
without tunes does not provide a satisfactory descrip- 
tion of the charged particle vs. J/^E* yield correlation, we 



employ the 



dN 



vs Nch distribution from |40| for the 



initial J/ 4* production. We now investigate the influence 
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FIG. 2: (Color online) Number of J/'i/ particles dependent on 
the charged particle multiplicity in proton-proton collisions. 
UrQMD calculations are compared to measurements at AL- 
ICE [U. We employed a rapidity cut of \y\ < 0.9 for the J/*s 
and ri < 1 for the charged particles to make our simulation 
comparable to the ALICE measurements. 



of the 'medium' created by analysis of the pp collisions in 
various charged particle multiplicity bins. Fig. [5]shows a 
comparison of our calculation to data measured at LHC. 
At low multiplicities one observes a good agreement to 
data. At higher charged particle multiplicities we ob- 
serve a slight suppression of the J/^* yield compared to 
the data. However, one should note that the magnitude 
of the input by [40] has been tuned to describe the pp 
data without additional late stage effects, meaning that a 
slight re-adjustment of the parameters might improve the 
description of the data. To make the J/ suppression in 
high multiplicity pp reactions more visible wc show the 
ratio Rpp of finally observed J/^'s over initially produced 
J/ 5's in Figl3] for each multiplicity bin at mid rapidity 

(|y|<i): 



|y|<i 



(6) 



At very high charged particle multiplicities in pp the 
J/4's might be suppressed by ^ 20% as compared to the 
initially produced J/4's at this multiplicity due to the 
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FIG. 3: (Color online) Number of J/^ particles dependent 
on the charged particle multiplicity. A rapidity cut of < 
1 is employed for all particles. The J/^P suppression in pp 
collisions at LHC and the Pb-Pb collisions at SPS show a 
very similar dependence on the charged particle multiplicity. 



phase has httle effect on the J/vt abundance. Let us 
also compare our results to the J/^ suppression in Pb- 
Pb collisions at SPS energies at similar charged particle 
densities (Fig. [31 squares). The similarity suggests that 
the J / 5* suppression may mainly depend on the charged 
particle multiplicity but not on the collision energy - at 
least if J/^* recombination can be neglected. 

To summarize, we have explored the potential for char- 
monium suppression in inelastic proton-proton collisions 
at y^Sivjv = 7 TeV as function of the charged par- 
ticle multiplicity. To this aim we have supplemented 
the UrQMD simulations with a phenomcnological pre- 
hadronic stage to allow for J/^ dissociation and melt- 
ing in the hot and dense medium. It turns out that a 
small amount of J/^ suppression (up to 20 — 30% at 
dNchI dy\\y\^i > 70) may be encountered even in pp in- 
teractions at LHC. 
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